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Abstract — In this paper, we consider a one-dimensional 
random geometric graph process with the inter-nodal gaps 
evolving according to an exponential AR(1) process, which 
may serve as a mobile wireless network model. The transition 
probability matrix and stationary distribution are derived for 
the Markov chains in terms of network connectivity and the 
number of components. We characterize an algorithm for the 
hitting time regarding disconnectivity. In addition, we also 
study topological properties for static snapshots. We obtain the 
degree distributions as well as asymptotic precise bounds and 
strong law of large numbers for connectivity threshold distance 
and the largest nearest neighbor distance amongst others. Both 
closed form results and limit theorems are provided. 

Keywords -random geometric graph; autoregressive process; 
component; connectivity; mobile network. 

I. Introduction 

Many randomly deployed networks, such as wireless 
sensor networks, are properly characterized by random geo- 
metric graphs (RGGs). Given a specified norm on the space 
under consideration, an RGG is usually obtained by placing 
a set of n vertices independently at random according to 
some spatial probability distribution and connecting two 
vertices by an edge if and only if their distance is less 
than a critical cutoff r. Topological properties of RGGs 
are comprehensively summarized in (TJ; also see (2) for a 
recent survey in the context of wireless networks. Although 
extensive simulations and empirical studies are performed 
in dynamical RGGs, analytical treatments of topological 
properties are merely done in static RGGs in the previous 
work. A recent paper is a remarkable exception, in 
which the authors conduct the first analytical research on 
the connectivity of mobile RGG in the torus [0, l) 2 . In this 
paper, we will also present analytical results and consider a 
one-dimensional exponential RGG process G(t, r, A) evolv- 
ing with time, where vertices are randomly placed along a 
semi-infinite line. One-dimensional exponential RGGs have 
been recently investigated by some authors |@)-|[6], which 
offer a significant variant from the familiar uniformly U[0, 1] 
distributed nodes, see e.g. 0, GJ-JU and references therein. 

'This paper is an extended version of a conference paper "Exponential 
Random Geometric Graph Process Models for Mobile Wireless Networks" 
pp. 56-61, presented in International Conference on Cyber-Enabled Dis- 
tributed Computing and Knowledge Discovery, 2009. 



In ifTOl . the distributions of distances between successive 
vertices rather than those of vertices themselves are exam- 
ined, and as it is stated in the same paper, this assumption is 
more natural since "sensors are usually thrown one by one 
along a trajectory of a vehicle^ We will then follow suit, and 
assume exponential distributions for inter-nodal distances 
of the graph process G(t, r, A). Every segment between 
two successive vertices is supposed to evolve following a 
stationary TEAR(l) process ifTTl with exponential marginal. 
This linear process has no zero-defect and thus surpasses 
the elementary AR(1) process involved in (6). We believe 
such a mobile scheme has broad potential applications due 
to the flexible double randomness mechanism (see Section 
II). Since the evolution of connectivity and the number 
of components in G(t, r, A) are both Markovian, we will 
address the transition probabilities and limiting distributions 
of these two processes Gt and G' t by employing Markov 
chain theory fl2l . |fl~3). It is worth noting that there are 
several Markov chains coupled in our model stemming 
from the first order autoregressive properties endowed in 
the evolution of inter-nodal distances. 

In addition to dynamical properties, we also establish 
static properties for fixed t. Vertices in G(t, r, A), for 
any given t, form nearly a Poisson point process (more 
precisely, a continuous time pure birth Markov process). 
The connectivity of a Poisson RGG is well-studied in the 
literature (see e.g. lfT4l - lfT7l ). especially in the context of 
ad hoc networks. We will investigate some topological 
properties basically along the lines of (4). In our opinion, 
the aforementioned simple idea in [ 1 1 reflects a conception 
of one step "memory" essentially. We show (in Theorem 4) 
that "1-step memory" + "growth" are not enough to produce 
power law distribution reminiscent of the architecture of 
Polya urn process, where typically infinite memory generates 
the power law lfl8l . 

Both finite and asymptotic analysis are given in this paper. 
We remark here that exact solutions are important since the 
asymptotic results can not be applied to real networks when 
not knowing the rate of convergence. 

The rest of this paper is organized as follows. Section II 
provides the definition of the exponential RGG process and 
some preliminaries. Section III deals with the evolutionary 
properties of G(t, r, A), including the transition probability, 



the stationary distribution and the hitting time for dis- 
connectivity. In Section IV, we present static topological 
properties of G(t, r, A) for fixed t. The degree distribution 
and strong laws of connectivity and the largest nearest 
neighbor distances are given among other things. In Section 
V, some concluding remarks and future research topics are 
discussed. 

II. Model and Preliminaries 

The RGG process G(t, r, A) is constructed as a discrete 
time process with n vertices deployed in one dimension 
on [0,oo). Let X\,--- , X* denote the vertices of the 
network at time t, for t > 0. Set Y* := Xf +1 - Xf, 
for I = 1, 2, • • • ,11 — 1 and Y" f := X\; see Fig.l for an 
illustration. We may envision time evolving upward along 
the i-axis and n vertices possibly growing along the x-axis. 

For < p < 1, we assume that {Y*} evolves following: 
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where the innovation sequences {eJ}t>o consist of i.i.d. 
nonnegative random variables. The behavior of this autore- 
gressive process {5j*}t>o is characterized by runs of rising 
values (with geometrically distributed run length) when 
choosing K* + e\, followed by a sharp fall when choosing 
e| without inclusion of the previous values. Furthermore, we 
assume that Yf, I = 0, 1, ■ • ■ , n — 1 are independent for any 
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Figure 1. One-dimensional exponential RGG process model 

III. Evolutionary Properties of G(t, r, A) 
A. Stationary Distribution of Gt 

Let us denote by Ct and T> t the events that G(t, r, A) is 
connected and disconnected at time f, respectively. Define 
Gt as a discrete time stochastic process describing connec- 
tivity of the graph process G(t, r, A). Therefore Ct — {Gt = 
"conneted"} and V t — {Gt = "disconneted"}. It is easy 
to see that Gt is a homogeneous Markov chain, assuming 
the cutoff r is independent of t. We abbreviate as usual the 
states as l="connected" (C) and 2="disconnected" (T>). Our 
main results in this section then read as follows: 

Theorem 1: Gt is a time-reversible, homogeneous finite 
Markov chain, with one step transition probability matrix 



P(n) 



Pn P12 
P21 P22 



where 



Pu 



1=1 v J 



(2) 



In particular, we set e\ :— (l—p)Zf, where Z\ ~ Exp(Xi) 
is an exponential random variable with mean A," 1 > 0. Let 
A := {Ao, Ai, • ■ ■ , A„_i}. In this case, as is shown in ifTTl . 
the above TEAR(l) process {^*}t>o would be a stationary 
sequence of marginally exponentially distributed random 
variables with parameter A;, assuming that the initial inter- 
nodal gaps Y[° are exponentially distributed with parameter 
Az. That means Y* 



l i 

correlation function of {Yf} is Corr(Y i t , Y, 



Exp(Xi). In this case, the auto 
'} is Corr(Y;*, Y { t+j ) = pi, being 
nonnegative. Reference |fl9l showed that ([T]) is stationary for 
each < p < 1 iff Yf is geometrically infinitely divisible. 
For further extension and discussion of ([T} we refer the 
reader to |20|. 

Remark 1: Vertices in snapshot of G(t, r, A) yield a 
counting process with inter-nodal distances having distri- 
bution Exp(Xi), while in standard exponential RGG, the 
corresponding distributions are relevant to n the total number 
of vertices (see Lemma 1 in flU) relying on the global 
information. 

Remark 2: Notice that the cutoff r = r(n, t) may depend 
on n and t. However, we restrict ourselves to fixed r in 
order to keep calculations clear though some results may be 
generalized without much effort. The popular assumption 
lim n _>oo r ( n ) — is not necessary here in virtue of 
unbounded support resulting from inter-nodal spacing. 
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(3) 

P12 = 1 - Pu and p 2 2 = 1 - P2i- 

Proposition 1: Gt has a unique stationary distribution 

7r(n) = (7Ti(n), 7T2(n)), where 

(I-P22) 2 

-P22) 2 +P2lPl2(2 -P22) fA . 
\2 



7ri(n) 



and 7T2(n) = 



(1-Pn) 2 



P22(l -Pll) 2 +Pl2P2l(2 -Pll)' 

Proposition 2: Suppose A; = A, for I — 0,1, ■ ■ ■ ,n — 1. 
Let P(oo) be the transition probability matrix of Gt as n 
tends to infinity, and 7r(oo) the (unique) stationary distribu- 
tion corresponding to P(oo). Then 7r(oo) = (0, 1) and 

lim n(n)P(n) = 7r(oo)P(oo). 

n — ^00 

Proposition 2 implies that we can swap the order of 
obtaining stationary distribution and taking limit w.r.t. n. 



Proof of Theorem 1: The probability density function of e\ 
can be shown to be given by fi(s) — Yz^e~ x ' s ^ 1 ~ p ^[ s >o]- 
Also, the conditional density function for Y* in the con- 
nected network is gy^civ) — r^r^xp l[o<y<r]. since the 
connectivity of network means Yf < r for all I. By 
independence property, we have p\\ = P(Ct+i\Ct) — 
I\i=i P(Xi +1 < f\Yi < r )- Our aim now turns to evaluate 
the probability P{Y l t+1 < r\Yf < r). Let V{ ~ Ber(p) 
independently, then the scheme ([TJ becomes 



if +1 = £ f + v^y/. 



(5) 



Let Y t t+1 denote Y{ +1 conditional on {Yf < r}. For a 
nonnegative random variable X with density function f(x), 
the Laplace-Stieltjes transform is defined by C(X)(s) = 
C(f)(s) = J °° f(x)e- sx dx. We have by ©, 

£(^ t+1 )( s ) = £( £ *)( S )-AW)( S ) 



l-p 



i-p di 



dy 



A, 



Xi + s{l-p) 



{(I-P) 



(s + Xi)(l-e-^' r ), 

where 5(y) is the Dirac-delta function. Inverting the above 
to get 

\ l (2-p)y 

C-\C{Y^))(y) = X ie -^l [v>0] + ^f-^- 

■ sinh v2(r^H> ] 

1 - e- A i' r 

• smh i 2(1 -p) 

Hence 

P(lf +1 < r\Y{ < r) = f C-\C{Yl +1 )){y)dy 

Jo 

(l-p)e- Xir (l - e~^) 
1 - e- x ' r 

(6) 

which gives (0. 

Let ^ A C [n - 1] := {1, 2, • ■ • , n - 1}. Denote the 
event E A := {Y* > r,Vl G A; Y,* < r,VZ G [n - 1]\A}, 



then we have 

p(c t+1 \E A ) = nm !+i <^>^) 

leA 

l£[n-l]\A 

= l[(l-p)(l-e-^) 



leA 



■ no- 

Je[n-l]\i4 

Here we used the expression P(Y l 



{l-p)e- Xir (l 



t+i 



-A,r 



< r\Y? 



> 



r) = (1 - p)(l - e^T^). Since P(£U) 
UieA^ Xir Uie[n-x]\A^ - e- Air ) and P(P t ) = 1 
Lir^C 1 - e ~ Air X © follows by noting that 



P21 =P{C t +l\V t ) = 



E 

07*AC[n-l] 



P(C t+1 \E A yP(E A )/P(V t ). 



Gt is time-reversible by standard results of Markov chains 

G2. ■ 

Proof of Proposition 1: Since Gt is an irreducible finite 
Markov chain, C and 2? are both positive recurrent. Also 
since they are both non-periodical, C and V are ergodic 
states. Set := min{/c : k > l,Gk = j,Go = i}, for 
i,j G {1, 2}. If the righthand side of the above definition is 
0, set Tij = oo. The first hitting probability is then given by 

= PWi = k\G =i). 
By a standard result from [13], an irreducible ergodic 
Markov chain has unique stationary distribution ir(n), and 
7Ti(n) is given by m(n) = l/Y,kLi ^Hi > for * = 1,2 in 
the present case. Thereby, follows easily from the facts 



In this case, 



/ff =Pii. =P2i^ 2 " 2 Pi2, for fc > 2; and /£> = p 22 , 
/ 2 2 =Pi2Pn 2 P2i, forfc>2. ■ 
Proof of Proposition 2: When A; = A, the righthand side 
of expression (O belongs to interval (0, 1). Hence pu tends 
to as n — > oo in view of (0. Since (1 — p)e~ Xr (l — 
e _ ^) + (l-e- Ar -(l-p)e- Ar (l-e - ^)) = l-e~ Ar < 
1, p 2 i tends to as n — > 00 by the binomial theorem 

and ©. Then we have P(oo) = ^ q j 

C is a transient state and V is an absorbing and positive 
recurrent state. By a standard result (see e.g. 0~3|), the 
stationary distribution corresponding to P(oo) exists and 
is unique. Direct calculation gives 7r(oo) = (0, 1). It is 
straightforward to verify that 7r(n) — + 7r(oo) as n tends to 
infinity. The theorem is thus concluded by exploiting the 
relation irP — it. ■ 

B. Transition Probability Matrix of G' t 

In this section, we show a refinement stochastic process 
G' t from Gt- To be precise, let {G' t = i} denote the event 
that G(t, r, A) has i components at time t, for 1 < i < n. 



Therefore, G' t is a homogeneous Markov chain with state 
space [n]. It's clear that {G' t = 1} = C t . 

Let the transition probabilities of G' t be p'^ := P(G' t+1 = 
j\G' t = i). Set A,BC[n — l] with \A\=i-l and \B\ = 
j — 1, 1 < i,j < n. Denote the event Ea ■= {Y( > r,VZ S 
A; F/ < r, V/ € [n. — 1]\A} and similarly for Pg. We obtain 
by the total probability formula, 

P«= E ^(^1^)-^)/^ = *), (7) 

A,BC[n-l] 

|A|=i=l,|B|=j-l 

for 1 < i,j < n. We have derived P(Ea) in the proof of 
Theorem 1, and P(G' t = i) = £ P(^a). To 

AC[n— l],|A|=i— 1 

evaluate (0, we still need the probability P(Eb\Ea), but it 
is also at hand already: 

P{E B \E A ) = [] P(Y[ t+1 > rlY* > r) 
leAnB 

■ J] P(Y t t+1 <r\Yt>r) 
l£A\B 

■ J] P(lf +1 >r|lf <r) 

ZG-B\,4 

• J] P(^ t+1 < r|F/ < r). 

Ze[ri-1]\AUB 

The second and fourth terms in the above expression have 
been obtained in the proof of Theorem 1, and clearly 



P(Y, 



6+1 



> r\Y, > r) 



1 



r),P(lf +1 > r\Y{ < r) = 1 - P(^ t+1 < r|lf < 



P(Y, t+1 < r|Y,' > 
t+i 



Now we arrive at the main result. 

Theorem 2: The transition probability matrix of G' t is 
P' = (Pij)nxn, which is given by ©. 

Of course, we have p' n — pn and 2j=2 Pij = Pi2- Since 
is an irreducible ergodic chain, it has a unique stationary 
distribution which may be deduced analogously as in Section 

m.A. 

C. Hitting Time for Disconnectivity 

Suppose Ct holds at time t, and we will consider 
the Markov chain Gt- Denote T := min{fc : fc > 
1, "Dt+k holds}, then T is the hitting time for disconnectivity. 
We may obtain the expectation of T using the transition 
probabilities derived in Section III.A by a routine approach 
iTPTl . In this section, we will instead depict an algorithm for 
getting the distribution of T directly. 

The event {T > fc} is equivalent to {Y t t+1 < r, Y,* +2 < 
r, • • ■ ,Y t t+k < r, V 1 < I < n - 1}. In view of (0, we can 
interpret the above as follows 



v t+2 



v t+k 



= e 



t+i-t 



j+k-i 



V, 



t+fc-1 t+fc-2 



+v, 



t+k- 



1 • • • Vf+'el 



V, 



t+k- 



1 • • • ViYi < r. 



Set U\ +1 := Vl +j e\ +j - 1 



Vi +j ■ ■ ■ VfYf, for 1 < j < k - 1 and XJ\ := VfYf. 
Therefore, conditioned on Y*, V*, ■ ■ ■ ,V l t+k ~ 1 , the prob- 
ability that the above k inequalities holds simultaneously is 
shown to be given by 

PhY l \{V l \---,V* +k - 1 }) = 

f Q - U > MeDdsj .... r;-^" 1 Me^)ds^\ (8) 

where /;(•) is given in the proof of Theorem 1. Denote the 
last i+ 1 integrals of dHJ by I^k-i, < i < k — 1. For i = 0, 



n,k 



rrt+k—l 



e 1 —p as = 1 — e l -*> 



1-p 



For i = l, 



h 



k-l 



r-Uf 



A, 



-e lLk<ie l 



t+k— 2 



= 1 — e l -p 

\,(r-Ut +k - 2 ) ^-»! +k - 2 ) 
\~p 6 1 [v i '+'- 1 =i] 

-11 -e i-p Je 1 -»1 [v h*-i =0] . 



In general, for < i < A; — 1, 



h,k-i 



r-Uf 



A, 



i,k-i+ide l 



.t+k-i-l 



io 1 -P 

We can proceed using this recursive formula by 
induction and integration by parts. Notice that 
P ; fc (^*, {VI, • • • , V^- 1 }) = 7 U from ®. 

Consequently, given Y* < r, the probability that Y* +1 < 
r, Y t t+2 < r, ■■• ,yj t+fc < r all are simultaneously true is 
seen to be given by 



Pf := 



A, 



1 - e- x ' 



E 



fc -vector g 
sting of fc l'a,fe— i 0' 



I' P*(y,0e- XlV dy. 
Jo 



Now we state our result as follows, whose proof is straight- 
forward at this stage. 

Theorem 3: Suppose the hitting time T of Gt is defined 
as above, then the distribution P(T < k) = 1 - n™^ 1 P i 
and it's expectation ET = Ilz^i 1 P k - 

In principle, by the truncation of fc, we may approximate 
ET arbitrarily close. 

IV. Snapshots of G(t, r, A) 

For fixed t, we denote by G(r, A) the static case which 
can be regarded as a snapshot of the dynamical process 
G(t, r, A). Also, we omit the superscript t typically, e.g. Yi, 
etc. 
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Figure 2. Probability that G(r, A) contains k components for different 
values of k 



A. Cluster Structure 

Let P n (C) denote the probability that G(r, A) is con- 
nected. We have the following result regarding connectivity. 
The proof is easy and hence omitted. 

Proposition 3: We have 



p n (c) = l[(i 



-Ajr 



)• 



8=1 



Moreover, suppose there exists A/ > such that A; < M, 
for all i, then P n (C) — » as n — > oo. 

Let ip n (k) denote the probability that G(r, A) consists of 
fc components and P™(fc) the probability that there are k 
components in G(r, A), each of which having size m (i.e. 
m vertices). 

Proposition 4: Suppose there exists M > such that 
A; < M, for all I. Then, for any fixed k, ip n (k) — > as 
n — » oo; and for any fixed k, m, P™(k) — > as n — > oo. 

Proof: Mimicking the proof of Theorem 3 and 4 in |4] 
yields the result. ■ 

In Figure 2, we plot ip n (k) as a function of n number of 
vertices for different k. We take A; = 1 for 1 < i < 10, 
and Aj = 2 for i > 10. Observe that the convergence to the 
asymptotic value is very fast. 

We may thus conclude that this static network is almost 
surely divided into an infinite number of finite clusters. This 
observation was first made in [15] by a different approach. 

B. Degree Distribution 

Let G(r, A) denote the graph G(r, A) when A = 
{A,-- - ,A}. 

Theorem 4: In the graph G(r, A), the degree distribution 
can be divided into three classes: the degree distribution of 

X\ and X n is Poi(Xr); and for fc + 1 < i < n—k, that of Xi 



T~}kefi' ^ or 2 < « < fc, the degree distribution 



is ,., j , ., 

of X, and X n+1 ^ is {e"»^ ©} 

Proof: Let {Yi}, {^/} be independent Exp(X). Denote 
the degree of vertex Xi as rfj. We get 



P(d n > k) 



P(di >k) = P{Y 1 



+ ifc < r) 



Xr f(Xr) k (Ar)^ 1 
fc! (fc + 1)! 



where we used an equivalent definition of gamma distribu- 
tion. Hence, 



P(d n = fc) = P(di = fc) = e 
Next, for 2 < i < fc, 

P{d n+X -i = k) = P{d t = k) 



fc! 



3=0 

Yt + --- + Y j+1 > r ) 
■P{Y{ + ■ ■ 
Yl + --- + Y>_ j+1 >r) 

.(As)'- 1 



< r, 



Ac: 



-Ax 1 



§7o'~ 0- - 1)! 

/•OO 

*/ r—x 

■J Xe Jk-J^ry- 

j Xe~ Xy <ly&x 

J r—x 



-2Ar 



(Ar) 
fc! 



E 

J'=0 



Finally, for k + 1 <i <n — fc 

fe 

P(d( - fc) 



n'-, < r, 



n + ■ • • + Y j+1 > r) 

■p(y; + ■ ■ 

Y{ + -.. + Yl_ j+1 > r) 
-2A,(2Ar) fc 
fc! 

which concludes the proof. ■ 

C. Strong Law Results 

Define the connectivity distance c n := inf{?' > : 
G(r, A) is connected}; and the largest nearest neighbor 
distance b n ■= maxi<i<„ naa.i<j<. n> j^i{\Xi - Xj\}. We 
derive asymptotic tight bounds for c„ and strong law of 
large numbers for b n , as n tends to infinity. 



Theorem 5: In the graph G(r, A), we have 



(i) 



lira sup - < 1 and liminf - > 1 a.s. 

n—too 2 Inn n->oo Inn 



(ii) 



v Xbn 1 

Urn — — = 1 a.s. 

n— >oc In n 



Proof: (i) Observe that P(c„ > x) < J2?=i e ~ A ' x = 
(n — l)e~ Xx invoking the Boole inequality. Let e > 0. Take 
x = x n = (2 + e) lnn/A in the above expression and sum 
in n, then we get 

oc oo 

J2 P (°n >*») <J2 n ~ (1+£) < °°- 
n— 1 n— 1 

By the Borel-Cantelli lemma, P(c„ > x i.o.) = 0. Hence, 
limsup„^ 00 < 1 almost surely. 

On the other hand, P(c n < y) = - e~ X,v ) 

: I - e-^)"- 1 . Take y = y n = (l- e) lnn/A, then 



-(1-6) 



n-1 



< oo. 



n=l 



We conclude that liminf 4^*- > 1 a.s. by using the 
Borel-Cantelli lemma again. 

(ii) By the independence of {Yi}, we obtain 

P(&„>x) = P( u^ 1 > x} n {F, > a:}} 

u{Yi > x} u {f„_i > x» 

n-l 

< ^ P(Fi_i > x) • P(K ( > as) 

i=2 

+p(Yi > x) + p(y n _x >x) 

= (n - 2)e- 2Aa; + 2e- Aa; . 
Take x = x n = (2 + e) Inn/ (2A), then we get 

oo oo 
n—1 n—1 

By the Borel-Cantelli lemma, lim supjj^^ < 1 almost 
surely. 

On the other hand, 

P(b n <y) - P{ ri'lz} {{Y l - l <y}yJ{Y l <y}} 
n{Yi < y} n {K„_! < y}) 

LfJ 



< ][[ p ( y 2i-i < y) • P(Y 2i < y) 



Arguing similarly as in (i), we can get liminf rwoo -^ B - > 1 
a.s.. This completes the proof. ■ 



V. Concluding Remarks 

This paper dealt with random geometric graphs in one 
dimension in which the vertex positions were evolving time. 
The critical assumption that this evolution was modeled by 
describing an evolution equation for the change in the inter- 
nodal spacing. We studied some dynamical as well as static 
properties and results were given for fixed n total number 
of vertices as well as n tending to infinity. 

It is worth pointing out that this paper is only a prelimi- 
nary step on the investigation of exponential RGG process 
models. The idea of considering spacings may be extended 
to high dimensions in the following way. Deploy X\ accord- 
ing to a probability density /, then place X 2 with the same 
probability density substituting the location of X\ for the 
coordinate origin, and so forth. We deem that the growing 
scheme would be an important alternative from the typical 
binomial and Poisson cases (TJ. Other interesting directions 
include examination of "multiple spacings", reinforcing 1- 
step memory to finite steps memory even to infinite one, 
which could be possible to result in power law degree dis- 
tributions. Since we only treat the limit regime for constant 
A;, how to deal with A; approaching infinity is our future 
research. 
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